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ABSTRACT 



Four decades have gone by since the discovery that the 
solar wind at 1 AU seems to exist in two relatively distinct 
states: slow and fast. There is still no universal agreement 
concerning the primary physical cause of this apparently 
bimodal distribution, even in its simplest manifestation at 
solar minimum. In this presentation we review and ex- 
tend a series of ideas that link the different states of so- 
lar wind to the varying superradial geometry of magnetic 
flux tubes in the extended corona. Past researchers have 
emphasized different aspects of this relationship, and we 
attempt to disentangle some of the seemingly contradic- 
tory results. We apply the hypothesis of Wang and Shee- 
ley (as well as Kovalenko) that Alfven wave fluxes at the 
base are the same for all flux tubes to a recent model of 
non-WKB Alfven wave reflection and turbulent heating, 
and we predict coronal heating rates as a function of flux 
tube geometry. We compare the feedback of these heat- 
ing rates on the locations of Parker-type critical points, 
and we discuss the ranges of parameters that yield a re- 
alistic bifurcation of wind solutions into fast and slow. 
Finally, we discuss the need for next-generation coron- 
agraph spectroscopy of the extended corona — especially 
measurements of the electron temperature above 1.5 solar 
radii — in order to confirm and refine these ideas. 

1. INTRODUCTION 

The intertwined nature of solar wind acceleration and the 
"coronal heating problem" has been known since Parker 
(1958) postulated a transonic flow solution made possi- 
ble only by the high gas pressure of the corona. Mariner 
2 confirmed the existence of a continuous supersonic so- 
lar wind in interplanetary space just a few years after 
Parker's initially controversial work (for a first-hand ac- 
count of the discovery, see Neugebauer 1997). Mariner 
also showed that the wind exists in two relatively distinct 
states: slow (300-500 km/s) and fast (600-800 km/s). 
The slow component was initially believed to be the 
"ambient" background state (e.g., Hundhausen 1972), but 
it was eventually realized that the fast component was 



in general more quiet and steady (Feldman et al. 1976; 
Axford 1977). The polar passes of Ulysses in the 1990s 
confirmed this revised paradigm (Gosling 1996; Marsden 
2001). 

In the 1970s and 1980s it became increasingly evident 
that even the most sophisticated solar wind models could 
not produce a fast wind without the deposition of heat 
or momentum in some form into the corona (e.g., Hartle 
and Sturrock 1968; Holzer and Leer 1980). It also was 
realized that the geometry of the flow — i.e., whether the 
magnetic flux tubes were radially expanding cones or su- 
perradially flaring trumpets — could have a significant im- 
pact on the mass flux and wind speed. This paper briefly 
surveys geometry-related explanations for the observed 
distribution of solar wind speeds, and presents a predic- 
tion for coronal heating rates in fast vs. slow flux tubes as 
a consistency check on these ideas. 



2. CORONAL SOURCE REGIONS 



Even after several decades of ever-improving in situ and 
remote-sensing observations, there is still no universal 
agreement concerning the full range of coronal sources 
of the solar wind. It is clear that strong connections exist 
between large coronal holes and the highest-speed wind 
streams (Wilcox 1968; Krieger et al. 1973; Noci 1973; 
Zirker 1977; see, however, Habbal & Woo 2001). The 
more chaotic slow wind, though, may come from a multi- 
plicity of source regions. Two regions that are frequently 
cited as sources of slow wind are: ( 1 ) boundaries between 
coronal holes and large streamers that undergo strong su- 
perradial expansion in the corona, and (2) narrow plasma 
sheets that extend out from the tops of streamer cusps 
(Wang et al. 2000; Strachan et al. 2002). However, dur- 
ing active phases of the solar cycle, there is evidence that 
slow wind also emanates from small coronal holes (e.g., 
Nolte et al. 1976; Neugebauer et al. 1998) and active 
regions (Hick et al. 1995; Liewer et al. 2004). During 
the rising phase of solar activity, there seems to be an 
abrupt (<; 6 month) change in the magnetic connectiv- 
ity between field lines in the in situ ecliptic plane and 
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the Sun (see Figure 5 of Luhmann et al. 2002). At min- 
imum, a large fraction of these field lines map into the 
high-latitude northern and southern polar hole/streamer 
boundaries, but at maximum nearly all the field lines map 
into low-latitude active regions and small coronal holes. 
The majority of the most recent transition time was not 
observed by SOHO because of its 4-month mission inter- 
ruption in 1998. 

The remainder of this paper is concerned mainly with the 
dichotomy between high-speed wind that emerges from 
the central regions of large coronal holes and low-speed 
wind that emerges from the hole/streamer boundary re- 
gions. More work is needed to apply the ideas presented 
below to other potential slow-wind source regions (e.g., 
strong-field flux tubes rooted in or near active regions; 
see Wang 1994). 

3. "GEOMETRY IS DESTINY?" 



There is a strong empirical relationship between the solar 
wind speed u measured in situ and the inferred lateral ex- 
pansion of magnetic flux tubes near the Sun. Levine et al. 
(1977) and Wang & Sheeley (1990) found that the asymp- 
totic wind speed is inversely correlated with the amount 
of transverse flux-tube expansion between the solar sur- 
face and a reference point in the mid-corona (see also 
Arge & Pizzo 2000; Poduval & Zhao 2004; Fujiki, these 
proceedings). As illustrated in Figures la and lb, the 
field lines in the central regions of coronal holes undergo 
a relatively slow and gradual rate of superradial expan- 
sion, but the more distorted field lines at the hole/streamer 
boundaries undergo more rapid expansion. It should be 
noted, though, that the eventual flux tube expansion (i.e., 
between the Sun and r — > oo) for polar coronal holes is 
likely to exceed that of the streamer edges, despite the 
opposite trend seen when the expansion factor / is mea- 
sured between R@ and a coronal source surface. 

Several potential explanations for the observed anticorre- 
lation between wind speed and flux-tube expansion have 
been proposed (see § 4). However, it is worthwhile to 
begin examining such a relationship from the standpoint 
of the the equation of momentum conservation along a 
solar wind flux tube: 

\ u I ar ar 

where, for a plasma dominated by protons and electrons, 
the effective one-fluid most-probable speeds are defined 
as cty/^ = ksiTp ll f±+T e )/m p and collisions and external 
sources of momentum are neglected. The function F(r) 
appearing on the right-hand side is defined as 

GM G 2 f r , 2 (2 1 df\ 

and /(r) is the dimensionless flux-tube expansion factor 
(which is proportional to B~ l r~ 2 measured along a flux 
tube; see also Kopp & Holzer 1976). 



Local extrema in F(r) satisfy the Parker (1958) critical 
point condition. Vasquez et al. (2003) found that only 
the global minimum in F(r) gives a sonic/critical point 
location that allows a consistent and continuous solution 
for u(f) over the full range of distances from the Sun to 
1 AU. For monotonically increasing expansion factors 
like those over the poles, F(r) tends to exhibit a single 
minimum in the low corona (r w 2 Rq). For streamer- 
like expansion factors that peak near the cusp, another 
minimum in F(r) appears at a height well above the cusp; 
this new point tends to be the global minimum. The lat- 
ter kind of flux tube — i.e., one that allows a more distant 
critical point radius — seems to correspond directly to the 
slow-speed wind measured in situ (see also Wang 1994; 
Bravo & Stewart 1997; Chen & Hu 2002). 

Figure lc shows the radial locations of minima in F(r) 
along individually mapped flux tubes that range from the 
pole to the edge of the streamer belt (see corresponding 
labels A — > D in the other panels). Eq. © was solved 
using the magnetic field model of Banaszkiewicz et al. 
(1998) and an isothermal corona (T p]l = T p ± = T e = 
1 .75 MK) for simplicity. The outer critical point appears 
only for field lines having latitudes at r — > oo less than 
about 23° above and below the equator. In more physi- 
cally realistic models that include radial and latitudinal 
temperature variations (e.g., Vasquez et al. 2003), the 
outermost minimum in F(r) is the global minimum, and 
thus as one moves from the centers of coronal holes to 
their edges, the critical point moves outwards abruptly 
from <J 2 Rq to 3-6 Rq at a latitude still rather far re- 
moved from the streamer cusp. 



4. HEATING ABOVE & BELOW THE CRITICAL 
POINT 

Why does the height of the critical point matter? Phys- 
ically, the critical or singular point (equivalent to the 
sonic point for a hydrodynamic pressure-driven wind) is 
the location where the subsonic (i.e., nearly hydrostatic) 
coronal atmosphere gives way to the kinetic-energy- 
dominated supersonic flow. 1 Whether the critical point 
lies above or below the regions where most of the energy 
deposition occurs is a key factor in determining the nature 
of the wind: 

1 . If substantial heating occurs in the subsonic corona, 
its primary impact is to "puff up" the scale height, 
drawing more particles into the accelerating wind 
and thus increasing the mass flux. Roughly, the in- 
crease in energy flux due to the heating can be bal- 
anced by the increase in mass flux, so that the even- 
tual kinetic energy per particle is relatively unaf- 
fected and the wind speed may not change (relative 



The idealized Parker critical point loses some of its mathematical 
importance when solving the time-dependent momentum equation (e.g., 
Suess 1982) or when including the effects of viscosity (Axford & New- 
man 1967). However, the critical transonic 'branch' remains the robust 
stable time-steady solution in nearly all models with varying levels of 
sophistication (e.g., Holzer & Leer 1997; Velli 2001). 
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Figure 1. (a) Idealized solar-minimum magnetic field configuration from the model of Banaszkiewicz et al. (1998). 
Selected field lines are labeled A^D in all 3 panels, (b) Superradial flux-tube expansion factors (normalized to f = 1 
at the solar surface) for 4 selected field lines, (c) Possible radii of the sonic/critical point computed from local minima in 
F(r), shown for a fine grid of field lines in the Banaszkiewicz et al. model (not all shown in [a-b]). Colatitudes of field 
lines (measured from the pole) at r = Rq and at infinity are plotted on the left and right vertical axes, respectively. 



to an unheated model). In some scenarios the mass 
flux increase can be stronger than the energy flux in- 
crease, and the asymptotic wind speed decreases. 

2. If substantial heating occurs in the supersonic 
corona, the subsonic temperature is unaffected and 
the mass flux is unchanged. The local increase in en- 
ergy flux has nowhere else to go but into the kinetic 
energy of the wind, and the flow speed increases. 

(Leer & Holzer 1980; Pneuman 1980; Leer et al. 1982). 
The above dichotomy is often modeled by changing the 
height at which the bulk of the energy is deposited, but 
it can also occur if the heating remains the same and the 
height of the critical point changes (as discussed in § 3). 

A natural link can be made between geometry-related 
changes in the flow topology and the heating-related 
changes in the wind. Wang & Sheeley (1991) proposed 
that the observed anticorrelation between u and / is a 
by-product of equal amounts of Alfven wave flux emit- 
ted at the bases of all flux tubes (see also earlier work 
by Kovalenko 1978, 1981). Near the Sun, the Alfven 
wave flux Fa is proportional to pVa{SV±) 2 . The den- 
sity dependence in the product of Alfven speed Va and 
the squared Alfven wave amplitude (SV±) 2 cancels al- 
most exactly with the linear factor of p in the wave flux, 
thus leaving Fa proportional mainly to the radial mag- 
netic field strength B. The ratio of Fa at the critical point 
to its value at the photosphere thus scales as the ratio of 
B at the critical point to its value at the photosphere. The 



latter ratio of field strengths is proportional to 1 //, where 
/ is the coronal expansion factor as defined by Wang and 
Sheeley. For equal wave fluxes at the photosphere for all 
regions, coronal holes (with low /) will thus have a larger 
flux of Alfven waves at and above the critical point com- 
pared to streamers (that have high /). 

To summarize, for streamers [coronal holes], more of the 
Alfvenic energy flux should be deposited below [above] 
the critical point. This effect is complementary to the 
change in height of the critical point discussed above; i.e., 
for streamers [holes] the critical point is high [low]. 



5. TURBULENT HEATING: FAST VS. SLOW 



One aspect of the Wang/Sheeley/Kovalenko hypothesis 
that needs further clarification is the link between an in- 
creased Alfven wave flux and increased coronal heating. 
Alfven waves can exert a dissipationless wave-pressure 
force that can accelerate the wind (e.g., Isenberg & Holl- 
weg 1982), but their ability to heat the plasma is less 
well understood. One idea that has received much re- 
cent attention is that low-frequency Alfven waves can 
be damped in the corona by undergoing a turbulent cas- 
cade from large to small scales. Here we present an em- 
pirically constrained model of Alfvenic turbulence and 
predict the contrast in extended heating that occurs be- 
tween a polar coronal hole flux tube and a near-equatorial 
streamer edge flux tube. 
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Cranmer & van Ballegooijen (2005) built a comprehen- 
sive model of MHD turbulence in a polar coronal hole 
flux tube. This model follows the radial evolution of the 
power spectrum of non-WKB Alfven waves (i.e., waves 
propagating both outwards and inwards along the flux 
tube) from the photosphere to 4 AU, and allows the tur- 
bulent energy injection rate (and thus the heating rate) 
to be derived as a function of height. The Alfven waves 
have their origin in the transverse shaking of strong-field 
(~1500 G) thin flux tubes in the photosphere. The bot- 
tom boundary condition on the wave power spectrum 
was derived from measurements of G-band bright point 
motions in the photosphere (e.g., Nisenson et al. 2003). 
Below the mid-chromosphere, where the bright-point 
flux tubes are isolated and thin, the linear wave proper- 
ties are computed using a generalized form of the kink- 
mode wave equations derived by Spruit (1981). Above 
the mid-chromosphere, where the flux tubes have merged 
into a more homogeneous network "funnel" (e.g., Tu et 
al. 2005), the wind-modified non-WKB transport equa- 
tions of Heinemann & Olbert (1980) are solved. 

Figure 2 shows a summary of results from the published 
coronal-hole model (corresponding to 9 = in Figure 
1) and the new streamer-edge model (corresponding to 
field line D). Below the transition region, the coronal- 
hole and streamer-edge models were assumed to be iden- 
tical except for the mass flux. Figure 2a plots the adopted 
wind speed for both flux tubes, as constrained by mass 
flux conservation, empirical electron densities (e.g., Sit- 
tler & Guhathakurta 1999) and the Banaszkiewicz et al. 
(1998) magnetic field model, modified by including the 
thin tubes and funnels at low heights. Also shown are the 
frequency-integrated Alfven wave amplitudes (6V±) as 
computed from the non-WKB wave transport equations, 
with turbulent damping as described below. 

Figure 2b gives the derived heating rate — expressed as 
the energy flux density Q per unit mass density p — for 
the two flux tube models. The slowly-varying ratio Q/p 
is plotted for convenience, because Q itself drops by 
more than 15 orders of magnitude between the transi- 
tion region (r s=a 1 .003 Rq) and 1 AU. In some ways, 
though, the plot is deceiving because it seems as if the 
streamer flux tube is heated more than the polar-hole 
model only below r 1.02 Rq. In fact, Q streame r > Qhoie 
everywhere below r w 1.4 Rq, but this is masked in the 
plot because ^streamer > Phoie- The heating rate is assumed 
to be equal to the energy cascade rate that has been de- 
rived for anisotropic MHD turbulence: 

Q - P (3) 

(Hossain et al. 1995; Matthaeus et al. 1999; Dmitruk et 
al. 2001, 2002), where £± is a transverse outer-scale cor- 
relation length and Elsasser (1950) variables are used to 
distinguish between outwardly propagating waves (ZJ) 
and inwardly propagating waves (Z+), with Z± = SV ± 
SB/y^Trp. The correlation length is assumed to expand 
with the transverse width of the flux tube (i.e., l\B = 
const; see Hollweg 1986), and its normalization is speci- 
fied by Cranmer & van Ballegooijen (2005). 
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Figure 2. In all panels, solid lines represent the coronal- 
hole flux tube (fast wind) and dotted lines represent the 
streamer-edge flux tube ( slow wind), (a) Outflow speed 
u and ID transverse Alfven wave amplitude (SV±) vs. 
height above the photosphere, (b) Heating rates per unit 
mass, (c) Estimated one-fluid temperatures correspond- 
ing to the heating rates in (b). UVCS-de rived proton tem- 
peratures and in-situ Helios one-fluid temperatures are 
shown for comparison (see text for references). 
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Figure 2c shows the result of integrating a one-fluid inter- 
nal energy equation (e.g., eq. 3 of Leer & Holzer 1980) 
to compute the mean temperatures T(r) that result from 
the heating rates discussed above. These results should 
be interpreted as preliminary because: (1) they are not 
the result of a self-consistent calculation of all fluid vari- 
ables, and (2) a simple choice was made for the electron 
heat flux (i.e., q e]] = qsn/10, where q$H is the classical 
Spitzer-Harm value) in order to split the difference be- 
tween the known strong conduction at low heights and 
the collisionless inhibition of q lle at large heights. How- 
ever, the overall trend in Figure 2c (i.e., the streamer 
being heated more at low heights, and less above the crit- 
ical point, than the coronal hole) is likely to remain valid 
when the above approximations are corrected. 

Several observations are also shown in Figure 2c, and the 
general agreement in the hole/streamer contrast lends cre- 
dence to the overall validity of the approach outlined in 
this work. UVCS/SOHO measured H I Lya resonance 
line profiles in the extended corona which provide a probe 
of proton velocity distributions. For off-limb observa- 
tions, the line of sight samples directions that are mainly 
perpendicular to the ^radial field lines, and the 1 /e line 
width V\ / e arises from two primary types of motion: 

2 = 2k^ + 

The two terms on the right represent random ther- 
mal motions and unresolved transverse wave motions. 
Using measured values of V\ i e and the modeled values of 
(SV±) for the waves, the above equation was solved for 
the plotted T p ±. The UVCS curve for a solar- minimum 
coronal hole (dark gray) comes from an amalgam of data 
from Kohl et al. (1997), Cranmer et al. (1999), Esser 
et al. (1999), Zangrilli et al. (1999), and Antonucci et 
al. (2000). The UVCS streamer data (light gray) came 
from Kohl et al. (1997) and Strachan et al. (2004). Also 
shown in Figure 2c are Helios data points which were 
computed by averaging together proton (Marsch et al. 
1982) and electron (Pilipp et al. 1990) measurements, 
and isotropizing from bi-Maxwellian fits; i.e., taking T = 
(T ll+ 2T ± )/3. 

When examining predictions for the plasma temperatures 
in fast vs. slow solar wind, it is worthwhile to compare 
with different, but potentially complementary ideas. Re- 
cently, Fisk (2003) and Schwadron & McComas (2003) 
discussed the origins of correlations between the eventual 
wind speed and observed properties of emerging loops 
in the low corona (see also the related footpoint diffu- 
sion model of Fisk & Schwadron 2001). Their predic- 
tion for there to be more basal coronal heating (and a 
higher mass flux) in the slow wind seems to be in accord 
with the results shown in Figure 2 (see also Matthaeus 
et al., these proceedings). There still seems to be a dis- 
connect, though, between theories of coronal heating via 
flux emergence and theories that invoke magnetic foot- 
point shaking (which in turn generates waves). The rel- 
ative contributions of these processes in various coronal 
regions needs to be quantified further. 



6. PROTON VS. ELECTRON HEATING 



The above analysis did not account explicitly for dif- 
ferences between heating the various particle species 
in the plasma. Even in a perfectly "collisionally cou- 
pled" plasma, there can be macroscopic dynamical con- 
sequences depending on how the energy is deposited into 
protons, electrons, and possibly heavy ions as well. 

Hansteen & Leer (1995) demonstrated these effects for 
a ID solar wind model: when all of the heat goes into 
electrons, there is substantially more downward conduc- 
tion compared to a proton-heated model. An electron- 
heated wind thus has a lower mid-corona temperature and 
a lower wind speed than a proton-heated wind. A 2D sim- 
ulation of streamers by Endeve et al. (2004) showed that 
the stability of closed-field regions is closely related to 
this kinetic partitioning of heat. When protons are heated 
strongly, the modeled streamers become unstable to the 
ejection of massive plasmoids; when the electrons are 
heated, the streamers are stable. Any sufficiently predic- 
tive model of fast and slow solar wind must take these 
effects into account (see also Cranmer & van Ballegooi- 
jen 2003). 



7. CONCLUSIONS AND FUTURE MISSIONS 



Our understanding of the dominant physics of solar wind 
acceleration has progressed rapidly in the SOHO era, 
and many of the insights embedded in the above anal- 
ysis would not have been possible without SOHO. In 
particular, the strong preferential heating and accelera- 
tion of heavy ions seen in coronal holes by UVCS has 
sharpened theoretical efforts to understand kinetic energy 
deposition in the collisionless extended corona (see re- 
views by Axford et al. 1999; Hollweg & Isenberg 2002; 
Cranmer 2002a; Marsch 2004). 

Despite these advances, the diagnostic capabilities of the 
SOHO instruments were limited and the most fundamen- 
tal questions have not yet been answered. If the kinetic 
properties of additional ions were to be measured in the 
extended corona (i.e., a wider sampling of charge/mass 
combinations) we could much better constrain the spe- 
cific kinds of waves that are present as well as the specific 
collisionless damping modes (Cranmer 2001, 2002b). 
Measuring the coronal electron temperature above 
~1.5 Rq (never done directly before) would allow us to 
determine the bulk-plasma heating rate in different solar 
wind structures, thus putting the firmest ever constraints 
on models of why the fast/slow wind is fast/slow. Mea- 
suring non-Maxwellian velocity distributions of electrons 
and positive ions would allow us to test specific models 
of MHD turbulence, cyclotron resonance, and velocity 
filtration. New capabilities such as these would be en- 
abled by greater photon sensitivity, an expanded wave- 
length range, and the use of measurements that hereto- 
fore have only been utilized in a testing capacity (e.g., 
Thomson-scattered H I Lya to obtain T e ). Spectroscopy 
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is key for the above measurements — especially in com- 
bination with coronagraph occultation — in order to mea- 
sure detailed plasma properties out into the wind's accel- 
eration region (see also Kohl et al., these proceedings). 
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